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Abstract 

We consider random flights of point particles inside ?2-dimensional channels of 
the form K*^ x B"~*^, where B"~*^ is a ball of radius r in dimension n- k. The particle 
velocities immediately after each collision with the boundary of the channel com- 
prise a Markov chain with a transition probabilities operator P that is determined 
by a choice of (billiard-like) random mechanical model of the particle-surface 
interaction at the "microscopic" scale. Our central concern is the relationship 
between the scattering properties encoded in P and the constant of diffusivity 
of a Brownian motion obtained by an appropriate limit of the random flight in 
the channel. Markov operators obtained in this way are natural (definition be- 
low), which means, in particular, that (1) the (at the surface) Maxwell-Boltzmann 
velocity distribution with a given surface temperature, when the surface model 
contains moving parts, or (2) the so-called Knudsen cosine law, when this model 
is purely geometric, is the stationary distribution of P. We show by a suitable 
generalization of a central limit theorem of Kipnis and Varadhan how the diffu- 
sivity is expressed in terms of the spectrum of P and compute, in the case of 2- 
dimensional channels, the exact values of the diffusivity for a class of parametric 
microscopic surface models of the above geometric type (2). 

1 Introduction and general definitions 

We consider mathematical models of particle-surface systems involving multiple, or iterated, 
(classical) scattering. The purpose of this section is to motivate the main question regarding 
these systems and to introduce a few definitions needed to informally explain our results. 
More detailed statements will be given in the course of the paper. 

1.1 An idealized experiment and the main question 

Figure fTTI depicts an ideal experiment in which a small amount of gas composed of point-like, 
non-interacting masses is injected into a (for simplicity of exposition 2 -dimensional) channel 
and the amount of outflowing gas per unit time is recorded. The graph on the right-hand side 
shows a typical exit flow curve. Possible gas transport characteristics that can be obtained 
from such an experiment are the mean value and higher moments of the molecular time of 
escape. 
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The central question then is: what can these time characteristics of the gas outflow tell us 
about the microscopic interaction (i.e., scattering properties) between gas molecules and the 
surface of the plates? 

For a more precise formulation of this question, we begin by describing the classical surface 
scattering operators that model the microscopic collisions of the point mass. We refer to the 
boundary of the channel region as the (wall) surface, irrespective of its actual dimension. 

Let H denote the upper-half plane, consisting of vectors {vi, V2) with positive second com- 
ponent. Elements of IH represent velocities of a point mass immediately after a collision with 
the surface. By identifying (i^i, - V2) and {vi, V2), we may regard pre-coUision velocities as also 
being in HI. A collision euentis then specified by a measurable map v eH^ rjp e CP(ll-ll), where 
T(ll-ll) indicates the space of probability measures on the upper-half plane. The measurability 
condition is understood as follows: For every essentially bounded Borel measurable function 
(p on HI, the function 

[P(p){v):^ \ (p{u)drjy{u) 
Ju 

is also measurable. We refer to P as the collision operator. This operator specifies the tran- 
sition probabilities of Markov chains with state space H giving the sequence of post-collision 
velocities from which the molecular random flight inside the channel can be obtained. 




Figure 1.1: Idealized experiment in which a small pulse of gas is injected into a 2-dimensional channel 
and the gas outflow is recorded. The graph on the right represents the rate at which gas 
escapes. From this function it is possible to derive mean exit time of escape t. The main 
problem is to relate easily measured properties of the gas outflow, such as t, to the micro- 
scopic scattering characteristics of the channel surface. 



If now T = T (L, r, s) denotes the expected exit time of the random flight, where s is the molec- 
ular root-mean square velocity, then a more restricted form of the general question is to un- 
derstand how T depends on P. (This expected exit time is easily measured in actual experi- 
ments involving gas diffusion using the rate of gas outflow as represented on the right hand 
side of Figure ll.ll see, for example, [22] for so-called TAP-experiments in chemical kinetics.) 

Although the analysis of P is generally simpler in dimension 2, an interesting complication 
arises here that is not present in the case of a 3-dimensional cylinder; namely, with respect 
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to the stationary distribution of velocities for natural collision operators (see the next subsec- 
tion), molecular displacement between collisions has infinite variance and standard central 
limit theorems for Markov chains do not apply. The same is true for the random flight in 
the region x [- r, r] between two parallel plates. In this regard, the random flight in a 3- 
dimensional cylindrical channel is simpler. (For an early study of two parallel plates case, see 
in f6l .) Infinite variance of the in-between collisions displacements requires a generalization 
of the central limit theorem of Kipnis and Varadhan in [20} that is proved in this paper. 
Now, from an appropriate central limit theorem we obtain for t the asymptotic expression 



(1.1) T{L,r,s) 
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for long channels in dimension 2, i.e., for large values of LI r, where D = D{r, s) is the diffu- 
sivity of a limit Brownian motion. Therefore, a more specific formulation of the problem is to 
understand how properties of the collision operator are reflected on D. A simple dimensional 
argument given in Subsection l2.2l shows that 

D(r,s) = 77, 

71 

where rj only depends on the scattering characteristics at the microscopic scale determined 
by P. The choice of constants will become clear shortly. 

The typical, but not the only type of operator we consider here is defined by a choice of 
microscopic contour of the channel wall surface, as suggested by Figure [LT] The main prob- 
lem then amounts to finding the functional dependence of 77 on geometric parameters of the 
surface microstructure. These parameters are scale invariant and are typically length ratios 
and angles. The presence of the logarithmic term in t is related to some surprising properties 
of D, as will be noted below, and for this reason we give somewhat greater prominence to the 
two-dimensional set-up in this paper. 



1.2 Natural collision operators and microstructures 

Let dV{v) denote the standard volume element on n- dimensional half- space H := H" and 
define the probability measure 

d^p{v)^27Ti^^^ ' {v,e„)expi^-^\v\^^dV{v). 

on HI. We refer to fip as the surface Maxwellian, or surface Maxwell-Boltzmann distribution, 
with parameter /3 and particle mass M. Here, (i;, e„> denotes the standard inner product (dot 
product) of v e H" and the unit normal vector, e„, to the boundary surface at the origin. We 
often denote this normal vector hy n- e„. It will be clear in context whether n refers to di- 
mension or to this normal vector. 

In physics textbooks, - IIkT, where T is absolute temperature and k is the Boltzmann 
constant. A simple integral evaluation shows that the mean squared post-collision ^peed with 
respect to /i^ is 

n C ■> n+\ 
Jh" pM 
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Another distribution of collision velocities that arises naturally is concentrated on an hemi- 
sphere S""" (5) := {f e H : I = s]; it is defined by 

r(^) 

d^iv)^ —{v,n)dV,^^{v) 

s"n— 

where dV,pi,{v) is the volume element on the hemisphere of radius s induced from the ambi- 
ent Euclidean space. In dimension 2, d^(v) - j^{v,n) dS{v), where dS indicates arclength 
element on S^{s). Equivalently, djiiO] - |cos0rf0, where 9 is the angle between v and the 
normal vector n. 

Definition 1 (Natural collision operators). The collision operator P will he called natural if 
one of the following holds: (a) fip is the unique stationary distribution for P, for some P; (b) 
the process defined by P does not change the particle speed and p. is a stationary probability 
measure for P for all s. If case (a) holds we say that the surface with associated operator P 
has temperature T - l/jf/3; in case (b) we say that P represents a random reflection. In addi- 
tion, we demand in both cases that P and its stationary probability satisfy the detailed balance 
condition (see equation (20.5) ofI23il]. We use v throughout the paper to indicate either pp or 




Figure 1.2: A periodic microstructure without moving parts. The cube containing a period of the mi- 
crostructure defines a cell; the point in a fc-dimensional torus T*' at which the particle enters 
a typical cell is assumed to be a uniform random variable. The particle enters with pre- 
collision velocity v and exits with post-collision velocity V. The corresponding P describes 
a random reflection, as defined in the text. 



The natural operators of particular interest to us are those defined by a choice of surface 
microscopic structure. We briefly described them here. (See [14], (TS), [16), and [17] for more 
details.) By surface (micro-) structure me mean that the channel wall's surface has a periodic 
relief composed of cells, each consisting of a mechanical system of moving masses or more 
simply a fixed geometric shape with no moving parts. (Natural collision operators specified 
by the former correspond to part (a) of the definition, while the latter correspond to part (b).) 
Moreover, the wall system, mechanical or purely geometric, is assumed to be at a "micro- 
scopic" length scale that, by definition, is incommensurate with that of the channel defined, 
say, by the radius of the ball factor B"~*^. At a collision event, the point particle enters a cell of 
the wall system, undergoes one or more deterministic collisions with it, transferring energy 
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between wall and particle in case (a), and leaves with a seemingly — from the perspective of 
the channel length scale — random velocity V . (See Figure [L2] ) Because of this assumption of 
incommensurability between the micro and macro scales, the relevant scattering properties 
specifying P are invariant under homotheties. Thus in dimension 2 we may, when convenient, 
assume that the width of each cell is 1. 

This incommensurability also dictates our assumption that the particle position on the en- 
trance of a cell at the beginning of a collision event is a uniformly distributed random variable. 
When the wall system has moving parts, the kinetic state of the cell at the moment the particle 
enters a cell is also drawn from a fixed probability distribution (a canonical Gibbs state at tem- 
perature r). Under these assumptions, V is an actual random variable and it can be shown 
(see (14)) that the associated collision operator P is natural according to Definition [1] The 
periodicity condition is not essential — all the basic facts discussed here hold, for example, for 
random structures, defined as probabilistic mixtures of periodic micro-structures. 

The operator P can be expressed as follows. Let / be, say, a continuous bounded function 
on H, and let 1/+ be the random velocity immediately after the collision of a particle with 
incoming velocity V- - v. Then 

{pf){v)^nf{v^W-^v\. 

For example, in the purely geometric case of random reflections in dimension 2, P is given by 

(P/)(0)= C f{^e{r^)dr 
Jo 

where ^e(r) is the angle V makes with the normal vector n and r e [0,1] is the position at 
which the particle enters a cell before collision. A similar integral over T*^ defines P in general 
dimension. 





mi 








mo 











Figure 1.3: An example of a wall system with moving parts. Mass niQ can move freely up and down, 
bouncing off elasticaUy against the fixed floor and an upper limit that is permeable to mi . 

Figure [L3l illustrates a wall system with micro -structure having moving parts. The periodic 
relief is assigned a mass mo and can move vertically and freely over a short range of distances 
[0, a\ from the fixed base, bouncing off elasticaUy at the lower and upper limits. The point 
particle of mass m\ enters the wall system with velocity y at a uniformly random location 
along [0,1]. Upon entrance, the velocity of the wall is assumed to be normally distributed 
with mean and a given variance a^, where mocr^ is proportional to the wall temperature; the 
particle then goes on to interact deterministically with the wall, leaving with random velocity 
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V . This is a very special example of wall system for which P is natural. For other examples and 
details about random billiards with microstructures omitted here see [14] . See also [18] for a 
detailed analysis of this and other examples in the context of stochastic processes in velocity 
space H having stationary measure \ip or /i. 

If we assume that the wall is static and has infinite mass, the system of the figure becomes 
purely geometric and the particle mass plays no role. In this case the operator P describes 
a random reflection. We refer to the subclass of natural operators derived from microscopic 
structures (either static or having moving parts, with arbitrary surface contours) as operators 
associated to surface microstructures. An interesting question suggested by (Tj is whether gen- 
eral natural operators are limits of operators associated to surface microstructures. 

Proposition 1. The operators for the classes of examples ofPisures il.2\ and\ 1.3\ are natural. 

Proof. See (141 [HI for proofs of these basic issues related to stationary measures and more 
examples. □ 

Occasionally, v will stand for either of the two measures fip or /i of Definition[T] Because 
P and V are assumed to satisfy the detailed balance condition, P is a self-adjoint operator 
of norm 1 on L^(lHI,v). A further assumption for the main results below is that P be quasi- 
compact, that is, the spectral radius of P restricted to the orthogonal complement of the con- 
stant function in (H, v) is strictly less than 1. Quasi-compactness for natural collision oper- 
ators for the types of systems illustrated in Figure [L2l is known to hold in a number of cases. 
The static version of the system of Figure [L2] (for the specific shape shown in the figure) has 
this property (see [16], [17]), and the operator for the one-dimensional version of the moving 
wall is known to be compact (see [14]); the case of the two-dimensional moving wall is still 
open. Further examples of shapes of systems of the static type having quasi-compact P will 
be provided later in this paper. 

Noting the two roles of P, as a self-adjoint operator on (H, v) and as a Markov transition 
probabilities operator, a useful characterization of powers of P is as follows. Let Vb, Vi, . . . be 
a stationary Markov chain with transitions P and initial distribution v and let 'P, <1> e (H, v) . 
Then it is not difficult to show that 

^vi',pfc<l)J>=Ev[^(V,-)0(T^-+fc)] 

for any i > 0, where Ey indicates expectation given that Vq is distributed according to v. 

2 Random flight and diffusivity 

2.1 BETWEEN-COLLISIONS DISPLACEMENTS AND TIMES 

The logarithmic term in Equation ll.ll is a special feature of the random billiard process in re- 
gions bounded by parallel plates in arbitrary dimensions (in particular 2-dimensional chan- 
nels bounded by a pair of parallel lines), and it is not present in the more typical cylindrical 
channel region IR*^ x EB"~*^ for k - 1, 2. Ultimately, this is due to the mean square dis- 
placements being infinite in the two-plates case and finite in the other cases, as will be seen 
later. This elementary but key observation is highlighted in the next proposition. 
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Let C :- C" :- x EB"~*^ denote the channel region. Of special interest are the low dimen- 
sional cases: n - k- 1, for « = 2,3 (two-dimensional channels and slabs in dimension 3) and 
A; = 1, n - i (cylindrical channels in dimension 3). Let represent the upper-half space 
consisting of vectors v e T^G, q e dC, such that (n, v) > 0, where n is the unit vector in per- 
pendicular to 56 and ■> is the inner product given by restriction of the standard dot product 
in R". If q' is the next collision point of the trajectory q+ f f , let Z(j;) denote the natural 
projection to the "horizontal" factor r'^ of the vector q' - qeR". We refer to Z(i') as the (hori- 
zontal) displacement vector for the given v. See Figure lOl The time of free flight between the 
coUisions at q and q' will be indicated by tj,(l'). 
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Figure 2.1: The between-coUisions displacement vector Z{v), where v is the post-collision velocity at 
a boundary point of the channel region. The time between two consecutive coUisions is 
denoted in this section Tj,(y). 
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Proposition 2. Letv he either of the two probability measures of Deflnition[^(denoted there pp 
and p). This is a probability measure on H = TqC (concentrated on a hemisphere in the case of 
p), for a given collision point q e 56. LetZa denote the product ofZ and the indicator function 
of the cone ¥\{a) :- {v e¥\: \Z[v)\ < ar] for a> 0. Also define for any unit vector u e R^ the 
orthogonal projection Zl^ :- (u,Za)- Then, ifn-k>2, 

E.[(Z'f|.limE.[(Z^f|.- 
Ifn - k-l, then the asymptotic expression 

Ey\{Z'jfj~4r^\na 

holds. The expected time of free flight tj, (here 'b' is for 'between collisions') is finite for both 
types of measures and n - k> 1. For the stationary measure p supported on the hemisphere of 
speed s, 

^^^^''^"s(«-fc) r(f) ■ 

Forpp, the corresponding expression is 



where M is particle mass. 
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For a sketch of the proof, see Section lSJ] 

Before considering random flights in C, we need to mention a technical point of geomet- 
ric interest that only arises in dimensions n > 3, having to do with whether the operator P 
actually gives rise to a well-defined process in 6. The issue is that in order for a fixed P to 
induce a scattering operator at each T^G, q e dC, we need to be able to identify the positive 
part (say, inward pointing) of this tangent space with the half-space H - {v e U." : (v, e„) > 0}. 
Such identification amounts to specifying an orthonormal flame on the tangent space at each 
boundary point, which provides the information of how the model microstructure is "aligned" 
with the channel wall. Introducing a frame field on 56 has, however, the effect that the scat- 
tering operator Pq at each boundary point q becomes a conjugate of P under an orthogo- 
nal transformation, rather than P itself. For all the cases in dimensions 2 and 3, namely 
In - 2,k - l),{n - 3,k - l),in - 3,k - 2), there is a natural frame (the parallel transported 
frame over the boundary of the channel) with respect to which Pq-P for all q. In the general 
case, TqQ has a canonical orthogonal decomposition into a "horizontal" part, naturally iden- 
tified with IR*^, and a vertical part that splits orthogonally into the normal direction Unq and 
the complement of dimension n-k- 1. We can now understand the issue as follows: Let v 
be a velocity vector of a particle that emerges from a collision qedQ and will collide again 
next at q' . In order that the Markov chain in velocity space be given by iterates of the same 
operator P, we need a field of orthonormal frames with respect to which v,atq, and its mirror 
reflection at q' have the same representation as vectors in IR". The components of these two 
vectors in Vj: and V"*-,, respectively, agree if we choose the canonical (parallel) frame, but on 
the subspaces themselves no such frame exists in general. With this in mind, and to avoid 
complicating the picture by introducing such frame fields as additional structure, we simply 
assume without further mention that P is V- isotropic, that is, it is invariant under conjugation 
by orthogonal linear maps that restrict to the identity on IR*^® [Re„. Notice that this assumption 
is vacuous in dimensions 2 and 3. 



2.2 DiFFUSIVITY, SPECTRUM AND MEAN EXIT TIME 

Let Xt, r > 0, be a piecewise linear path in the channel region 6 describing a random flight 
governed by a natural collision operator P. Recall that P is a self-adjoint operator on L^(1H1, v). 
We assume throughout that P is quasi-compact. (This is one of the conditions needed for 
Theorem[2l below.) 

We wish to consider a diffusion process in K*^ obtained by an appropriate scaling limit of the 
projection of Xt to the K*^ factor of C. Let this projection be denoted Xt, and assume Xq - 0. 
The sequence of post-collision velocities of Xt is a stationary Markov chain Vo,Vi,..., with 
initial distribution v. The displacement vectors, previously defined, are random variables 

Zo,Zi, Thus Xf , at collision times, are sums of the Z, . The a -scaled random flight is defined 

as follows. Let hia), a > 0, be 

{a for n - fc > 2 

a/ log a for n - k- 1. 

Define the scaled channel systemwith scale parameter a > to be the channel system with ra- 
dius rla and root-mean-square velocity h{a)s,ns- (If v = /i, Sms is the constant speed through- 
out the process and if v = /i|3, = {n+\)l pM.) The random flight paths and their projection 
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are defined as the paths for the a-system. We denote them by Xa,t and Xa,t, respectively. The 
a-scaled free displacement with post-collision velocity vector v is a~^Z{v), and the displace- 
ments associated to the Vj are a~^Zj. For any t > 0, over a time interval [0,t], the number of 
collisions of X t with 56 will written N-^ . For the a-scaled system this number is Na,T Nah{a)T ■ 
For the cases in which n-k-\, when Ey \\Z\^\ is infinite, it will be necessary to also consider 
the a-truncation Za of Z introduced above in Proposition[2] The a-scaled a-truncation of Z 
and Zj wiU be written a~^Za and a'^Z^j- Finally, we will like to follow the projected random 
flight along an axis set by a unit vector u e K^'. Thus we define Z", X", Z^', etc., to be the 
orthogonal projections of Z, Xr, Za, etc., on Uu. 

Theorem [1] below gives conditions under which Xa,t converges to Brownian motion for 
large a. In this subsection, we wish to focus on the variance (or diffusivity) of the limit Brow- 
nian motion, and provide an interpretation of this constant in a way that does not make use 
of the physically somewhat artificial scaling just introduced. The precise conditions for the 
diffusion limit to exist when Z has infinite variance relative to v (i.e., for « - A; = 1) are not yet 
fully clear; for Theorem[T]we make the following additional assumption that will be verified 
in the examples discussed later. 

Assumption 1. Define forY> 1 the setW'^ia} :- {v eH: \Z"iv)\ < arllog^ a} and let Z" . be 
the product ofZJ by the indicator function of W[a]. Then, for any t > 0, 



lim Ev 



I Nr,. 



-T Z" 

V " =0 



exists for all unit vectors u e [ 



This assumption, which will only be needed when n - k- 1, will be explained later in the 
context of Theorems [2] and |3l It is not needed for n - k>2, when Z has finite variance and 
Theorem[T]is for the most part a consequence of well known limit theorems, in particular the 
central limit theorem for reversible Markov chains as formulated by Kipnis and Varadhan in 

Ho]. 

Theorem 1 (Diffusion limit). LetP be quasi-compact and, ifn-k- 1, suppose that Assumption 
\l\holds. Then the a-scaled projected random path X^^, for a unit vector u e U'^, converges 
weakly as a ^ oo to a Brownian motion in Uu with diffusion constant D" further specified 
below. In particular, X^^ converges in distribution for each t>0 to a normal random variable 
in K.U with mean and variance tD". The following statements concerning Ti" also hold: 

1. For the purpose of having a baseline valueforD", suppose thatP maps probability mea- 
sures to V. In other words, let the velocity process on H be i.i.d. with probability measure 
V. Let n-k>2. Then, denoting the constant Dg in this special case, 



when v-pp and 



V2n{n+l) [n-kV-l 



g^g^A^m .rs 



0?r(^) («-A;)2-l 
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when V - 11. Recall that Sms is the root-mean-square velocity for fj.^, and that fi is con- 
centrated on the hemisphere of radius s in H. Being independent ofu, we denote these 
values by Dq . 

For n-k- \, the baseline diffusivities are 

_ 4 
y/2n{n+l) 

when v-pp and 



when v-p, supported on the hemisphere of radius s. 

3. The diffusion constant for a general P can now be written asT)" - ri{u)'Do, where rjiu) 
has the following expression in terms of the spectrum of P . First consider the case n-k> 2 
and define a probability measure on the spectrum by 

n"(dA) := \\Z"r^(Z",n{dViZ"), 

where U is projection-valued spectral measure associated to P and the inner product and 
norm are those ofL^ (H, v) . Then 

fi 1 + A „ 

(2.1) T][u)^ - — -n"WA). 

J-l I- A 

Notice that rjiu) is quadratic in u. Now suppose n - k - I and define for each a the 
probability measure 

nlidX) := \\z^\\-\z'J,nidA)z^) 

on the spectrum and the function r\a[u) - f\ j^YV^idX) . Then the limit lima^caila(u) 
exists and defines a quadratic function r]{u) ofu. 

Proof. The limit theorems in probability theory we require are standard in the case n - k>2 
and will be proved later for n - = 1 in Theorems |2]and|3] Here we only indicate how Do and 
the expression for 77(m) in terms of the spectrum of P are obtained. 

Recall that Na^t - Nah(a)t is the number of collisions of the a-scaled random flight with the 
boundary of the (a-scaled) channel region during time interval [0, f] . Let tj denote the time 
duration of the step with (non-scaled) displacement Zj. Then, as 

(To + --- + TiVj._i)/Afr< TINt<[tq + --- + tnj)INt 

for any T > 0, we can apply Birkhoff's ergodic theorem to obtain Ev[tj,] - limj-^oo T/Afr- 
where we have used a previous notation tj, for the random time between consecutive col- 
lisions. 

Although not necessary in this case, we use here the truncated displacement Za, so that 
the derivation of the spectral formula will also apply to the infinite variance case to be dis- 
cussed later (under the more stringent conditions needed in that case). It wiU be shown later 
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in Proposition[7]that, for any f > 0, 



(2.2) 

In the i.i.d. case, this gives 



D"= lim 4-Ev 



INa. 



D" = lim Ey Z" - hm Ey Z" ^ = hm 



h{a) 



jh{a)t 



h{a) 



Ev[Tb] 



We can now invoke Proposition|2]to obtain the values claimed for D" in the i.i.d. case. 

Next we obtain the spectral formula for t]{u), beginning from expression 12. 2 1 This expres- 
sion holds without further assumptions in the finite variance case, and it follows from As- 
sumption [T] when n - fc = 1 as will be shown later in Proposition |7l Because P has positive 
spectral gap and Z^ has zero mean, the measure n has compact support in the interval (-1,1). 
In particular, 1 - A is bounded away from zero on the support of (Z^ , n(dA)Z^). Now observe 
that, for i> i, 



\\z^fn"(dX]. 



a,j a,i\ 

With this in mind, we obtain for a fixed N after some algebraic manipulation. 



Ev 



(N-l 

\j=o 



N-lj-l 

N+2 ^ ^A-'- 

j=i i=o 



Iz^'fn^WA)^ 



1-hA 



■IN+om] wz'Jfn'iidA). 



The expectation on the right-hand side of limit lz!2l can be written as 



1 oo 



/iV-1 



Na,t+om 

ah[a)t 



^WZ'^fn^m). 



Keeping in mind the relationship between the expectation of Na,t and EylTt,] observed above 
in the derivation of the i.i.d. case, we have 



lim ^-Ey 



INa. 



E zlj 



= lim— — ^ — -n^WA) = Do — — n"(dA). 
fl-oo a Ev[Tb] j-il-A " j-il-A 



This proves the claimed form of the diffusion constant. The necessary central limit theorem 
and weak invariance principle required to prove convergence to Brownian motion (in the case 
n-k- 1) will be shown later. □ 

We remark now on a simple interpretation of the diffusivity D in the context of the idealized 
experiment described earlier that allows us to obtain D without recourse to the somewhat 
physically artificial a-scaling. In order to keep the discussion simple, only the case fc = 1 is 
considered, although the main idea can be generalized in obvious ways. 

Consider a channel C(L) - [-L,L] x EB"~^ (r) of length 2L and recall that t(L, r, s) is the mean 
exit time from C(L), introduced in Subsection ll.il where 5 is root-mean-square speed. The 
following elementary dimensional properties are easily derived: 
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(i) t(L, r, s] - T(aL, ar, as) 



(ii) t(L, r, s) - aT{L, r, as) 

(iii) T{aL,r,s) - ah{a)T{L,r I a,h{a)s) 

where the third property is a consequence of the first two. It also follows from (i) and (ii) 
that the function F{LI r) :- [si r)T{L, r, s) is independent of s, dimensionless (that is, devoid of 
physical units), and scale invariant. 

We are interested in the asymptotic behavior of t(L, r, s) as L grows to infinity. Since the 
mean exit time from the interval [- r, r] for Brownian motion with diffusivity D starting at is 

IT) we expect, given Theorem[T]and the above properties of the mean exit time. 



(2.3) T(L,r,s)' 



^ \fn-k>2 

'D\n{Llr) ^^n-k^\. 



Notice, in particular, the expected relation 

'D^C[P)rs, 



where C[P) - limo^oo ah{a)IP{a), being independent of L, r, s, is a characteristic number of 
the scattering process at a microscopic scale. This asymptotic expression is indeed true, and 
it is a consequence of the following proposition, which will be proved later. 

Proposition 3. Let L > and T be the function on the space of continuous paths j : [0, oo) K 
defined bylij) :- mf{t > : l7(f)| > L], where the infimum of the empty set is taken to be -oo. 
LetE^ denote expectation with respect to the law of the process t ^ Xa,t, conditioned to start at 
andEg, similarly defined, for the Brownian motion with diffusion constant D. Then 

lim E^[T] =E^[T] = 

a— oo 

and the asymptotic expression \2.3\ holds. 



2.3 Generalized Maxwell-Smoluchowski models 

Before examining diffusivity in examples of collision operators derived from geometric mi- 
crostructures, we briefly mention an interesting generalization of a very classical example 
widely used in kinetic theory of gases, known as the Maxwell-Smoluchowski collision model. 
The collision operator for this model is P^s - oiQp + (.1- a)I, where a e [0, 1] is a constant, I is 
the identity operator, and Qp is the projection into the subspace of constant functions: 

[Qn(p)iv):^ I (p[u)dp.a{u). 

JM 

The interpretation is that, upon collision with the surface, a particle scatters according to p^ 
(diffuse scattering; see definition of pp at the beginning of Subsection II. 2> with probability 
a, and reflects specularly with probability I - a. Clearly, is a natural collision operator 
according to Deflnition[T] A similar definition can be made for random reflection operators. 
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with /i in place of fip. In either case, the diffusivity is easily shown to be D = jz^Dq, where 
p- 1 - a is the probability of specular reflection and {1, p} is the spectrum of Pms- This simple 
model is very useful in providing a rough interpretation of the typical 77 obtained from micro- 
structures. 

A generalization is as follows. Let Q be a collision operator, not necessarily natural, such 
that for V eH the measure 

A^Q(v,A):^{QTiA)(v) 

is absolutely continuous with respect to the standard area measure, where ^4 is a measurable 
set and denotes its indicator function. We say that operators such as Q are diffuse. Let 
a : H X H — ► [0, 1] be a measurable function. Now form the Metropolis-Hastings kernel with 
proposal Q and acceptance probability a. The corresponding operator is 

[P(p] [v] - q){v) + I a{v,u)[(p{u) - (p(v)] Q{v,du]. 
Jh 

The interpretation is that, for an incoming velocity v, a candidate outgoing velocity V is cho- 
sen according to Q{v,-); then, with probability a{v,V), V is accepted as the post-collision 
velocity, and with probability I - a{v,V) the post- collision velocity is taken to be v itself. 

The standard Metropolis-Hastings method provides an explicit a for a given Q such that the 
resulting P is natural with respect to ^p. See (26| for the general construction. This provides a 
large class of examples of natural collision operators satisfying the following definition. 

Definition 2 (Generalized Maxwell- Smoluchowski models). We say that the collision operator 
P defines a generalized Maxwell-Smoluchowski model if it is a natural operator for a given 
IJ,p and is of the Metropolis-Hastings type with a diffuse proposal operator Q and acceptance 
function a. 

Two families of examples of generalized Maxwell-Smoluchowski models are given in the 
next section. In the family where a flat transition of (scale-free) length h is added between 
semicircles, the proposal Q is simply the operator associated to the semicircle geometry and 
the acceptance function is the constant 1 - ^. In the family where the semicircle is split and a 
flat floor of length h is added, the acceptance function is more complicated as it depends on 
the incoming pre-coUision angle. 

Regarding these operators and their diffusivities, we only indicate the following very gen- 
eral comparison result. Let Pi and P2 be natural collision operators with the same stationary 
distribution /i. We say that Pi dominates Pz off the diagonal ilPi{v,A\{v})> P2[v,A\{ v}) for 
|i-almost all v. The next proposition holds with little modification to what is given in f26]. 

Proposition 4. Let P\ and P2 be natural collision operators and let Di and D2 be the asso- 
ciated diffusivities for the limit Brownian motion process in the scaled channel system. If Pi 
dominates P2 off the diagonal then T)\ <D2. 

In other words, an operator that is more dispersing in the sense just defined has slower dif- 
fusion. This result is illustrated by the example of Figure l2.3l of the next section featuring a flat 
transition between semicircles. The formula forVih) (next section) shows that the diffusivity 
increases with h, which is what is expected from the proposition. A similar comparison of 
operators does not hold for the example shown below where a flat floor is added. 
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2.4 Examples of diffusivity for geometric microstructures 



We limit our attention in this subsection to operators associated to static microstructures in 
dimension 2. Therefore, the main question of interest in how the shape of the surface con- 
tour influences the signature parameter rj-VIVo of the particle-surface interaction. For the 
examples given here, rj can be obtained exactly. It will be noticed that the examples are vari- 
ations on a theme: they are built out of arcs of circle and straight lines and correspond to 
focusing billiards. 



Figure 2.2: A geometric microstructure consisting of semicircular arcs. 



(2.4) 



For the focusing semicircle structure of Figure lz!2l it will be shown that 

l-|log3 



l+ilog3' 



This value is reminiscent of the limit variance seen in central limit theorems for the stadium 
billiard of deterministic billiard dynamics (see [4]). 

A simple modification of the semicircles contour is shown in Figure 12.31 It consists of 
semicircles as in the first example separated by flat sections. We introduce the parameter 
h = II {I + 2r) e [0,1), which gives the proportion of the top line occupied by the flat part. 




Figure 2.3: Circles and flats. For the non-periodic shape on the right, we assume that the fraction of 
length comprising the flat part on top is well defined and equal to h = I! (l + 2r), where I and 
r are described on the left figure. The parameter 77 is the same in both cases. 



The diffusivity for the shapes of Figure [Z3l as a function of h, is given by 

(2.5) D(/i) = Do-/— - 

1- h 

where 77 is the signature diffusivity parameter of the example of Figure 12.21 Clearly, for the 
h-0 limiting case, the microscopic cell is simply the semicircle and D{0) - D. At the other 
end, as h approaches 1 the diffusivity increases without bound. Of course, for a completely flat 
surface, the transport ceases to be a diffusion at aU and becomes (a much faster) deterministic 
motion. 
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The next example refers to the surface of Pigure l2.4l In this case, the parameter h measures 
the length of the middle wall relative to the period length of the contour. 




Figure 2.4: Semicircles with middle wall. Define the scale free parameter h = II [2r). 



For the middle wall of relative height h<\l2 and h- 1/2, the values are, respectively. 



(2.6) 



l-7log3 



D(/z) = Ho T— ^, 11(1/2) = 2)0 ^ 



l + 4log3 



l + ilog3' 



l-ilogS' 



Observe, in particular, that the diffusivity does not change, and has the same 77 as the example 
of Figure 12.21 until the middle walls reach the top of the cell. At that point the diffusivity 
changes discontinuously to 2) (1/2). 

A related phenomenon is seen in the next family of examples, shown in Figure 12.51 It is 
obtained from the first example by adding a flat floor of relative length h- II {I + 2r) e [0, 1). It 
will be shown for this parametric family that 



(2.7) 

where 



l-(h' 



l + 3hl-h 3+h 
■log: 



1 + h 



1-h 




Figure 2.5: A family of geometric microstructures with parameter h = ll(l + 2r). 



At the h - limit we naturally have the same rj as for the first example. What happens 
when h approaches 1 is perhaps more surprising. In this case Vih] approaches the baseline 
value Do- Recall that this is the diffusivity of the process where at each collision event the 
particle reflects, independent of the pre-collision angle, according to the stationary measure 
/i. From the perspective of a single collision event, collisions are nearly mirror-like; on the 
other hand, from a multiple scattering perspective the collision process reaches equilibrium 
instantaneously making the surface ideally rough in a sense. This peculiar phenomenon and 
the discontinuity in D seen in the previous example are due to the fact that the diffusivity 
is determined only by collisions which occur at angles nearly parallel to the channel walls, a 
result made explicit in Proposition[7] An allusion to this property is found in (6j . 
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The computation of T) for these examples will be given in SectionH) 




Figure 2.6: The family of geometric microstructures of Figure [231 includes the above when I is allowed 
to be negative. Let Z = s - 2r, where r is the radius of the arcs of circle and s is the width of a 
period cell and 5 is the length of the opening. For l>-r, the diffusivity D(h), parametrized 
by h = lHl + 2r],is still given bvEQuation l2.7l 



3 The main limit theorems in the infinite variance case 

Let P be a natural collision operator v^^ith stationary measure v. Define for each y e IHl the 
measure P{v,A) = (PIa) (v) for j4 c IHl measurable. Let O = and denote by 3^ the product 
Borel (T-algebra. Define a measure P on cylinder sets as 

P{{a)en:a)oeAo,...,co„eA„})= I j ■■• / P[Li)n-i,da)n)--- Pim,da)i)fi{da>o), 

JAo JAi Ja„ 

for j4o, . . . , A,i c (H measurable and extend it to IT. The coordinate projections V/ : Q ^ HI given 
by Viiw) - (jLii, for / = 0, 1,... , thought of as random variables on the probability space (n,IT,P), 
constitute a Markov chain. Because v is stationary for P, it follows that the Markov chain is 
stationary. 

Being a natural operator, P together with v satisfy the detailed balance condition, which 
implies it is self-adjoint on L^(ll-ll,v). We also assume from now on that P is quasi-compact, 
and so has a spectral gap. These assumptions imply that the chain is p-mbdng. That is, 

sup {corr(X, Y):XeL^ (iTo^) ,YeL^ , A; > if = O (p"] , 

for some p such that < p < 1, where the notations are as follows: corr(X, Y) is the correlation 
between X and Y, (Jg) is the space of square integrable functions measurable with respect 
to the (T-algebra 7^ generated by Vo,...,Vic and 1^(3^^^) is the space of square integrable 
functions for the cr-algebra generated by Vj, for j >k + n. 

Consider the a-scaled channel system, as defined earlier, with channel radius scaled down 
by a and each post- collision velocity scaled up by the previously defined h{a). Let Xa,t be the 
position at time t along the horizontal axis of the particle for the a-scaled random flight start- 
ing at Xafl - 0; let Na,t be the number of collisions with the walls for the scaled system during 
the time interval [0, t\. Notice that Na,t - Nah(a)ty where Ns is, by definition, the number of 
collisions during [0, s] for the non- scaled system. 

In what follows, h[a) - a/ log a and Z will be any real measurable function onH slowly vary- 
ingat infinity, in the sense that Ey \Z^\z\^a\ ~ Clog a as a ^ oo, where C > is a constant and 
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a„/bn means a„lbn ^ 1 as n ^ oo. The various notations used before for the inter- collision 
displacement function will be used for this general Z. Thus, for example, Zj - ZiVj). How- 
ever, the ^-truncation Z„ of Z will be understood more generally as follows: If I{a) is an inter- 
val whose endpoints are functions of the scaling parameter, we write z'^"'' :- Z'i[Z£i{a)}y and 
if there is no ambiguity about which interval is being assumed we write Za '■- Z^'"\ Com- 
bining notations, Zaj - ZaiVj), for the non-scaled observable Z, while the corresponding a- 
scaled quantity (for the a-scaled system with radius r/a and root-mean-square speed hia)s) 
is a~^Za,j. Unless explicitly stated otherwise, Za is associated to the interval I{a) - [-a, a]. 
The probability measure UaidA) on the spectrum on P is defined as before for a general Z: 
YlaldX) := \\Za\\~^ {Za,Tl{dX}Za) ■ While we make no claims on the existence of a weak limit 
of the measures Ha we note that there exists a subsequence of Yla that converges weakly to a 
probability measure which we will call Do. 

If a sequence of random variables converges in distribution to a normal random variable 
with mean and variance cr^, we say for short that the sequence converges to 3\f (0, cr^) . 

Theorem 2 (Central limit theorem). Suppose P is a quasi- compact natural operator and that 
Assumptior^holds. Then, for any t>0, the sumY.^^'^ ^ converges fo3\f(0, fD) asa^oo, 

where 



•D= lim 



a— oo a t 



INa. 



\ 7=0 



Theorem 3 (Weak invariance principle). Under the same assumptions as in Theorem[l\ let 
Xa,t E IK*^ be the particle at time t in the a-scaled system with radiusri a and root-mean-square 
velocity h{a)s. Then Xa,t converges weakly to Bt, a Brownian motion with dijfusivity given by 
the quadratic form D". 

We note that it is still possible, without Assumption[T] to prove a central limit theorem and 
weak invariance principle for the inter-collision displacements — only slight modifications of 
the statements and proofs are needed. However, Assumption [1] allows us to express the vari- 
ance of the limit distribution in the central limit theorem in terms of the observable Z and 
the operator P. So while more general statements can be made and proven, we choose when 
possible to emphasize the connection between macroscopic data — the limit variance in the 
central limit theorem — and microscopic data encoded in the operator P and its spectrum. 

We also note that Assumption[T]may be reduced to a statement about covariances. Observe 
that 

E a-^Za,j 
\ i=a 

As already seen in the proof of Theorem [1] the limit of the first summand exists in general. 
Thus the assumption may be restated as requiring the existence of the limit of the second 
summand of covariances. In the case of strictly stationary p-mixing sequences with finite 
variance and a relatively light (less than exponential) condition on the rate of mixing, the 
scaled limit of variances is known to always exist (see for example 171). As far as we are aware, 
the corresponding result in our setting — whether such a limit always exists for exponentially 
fast p-mixing stationary sequences with infinite variance — has not been addressed. 



\2 




Na.t-\ 








= E 


i=o 


+ E 


O- E ^a,iZa,i 
0<\i-j\<Na,t 
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The remainder of the paper is organized as follows. Section ISTTI oudines the proof of the 
central limit theorem as a sequence of technical lemmas, leaving the proofs of the lemmas 
for Section [5^ In Section |3^ we prove the weak invariance principle and Proposition[3]on 
convergence of mean exit times. The computation of diffusivity for the examples of Section 
I2.4l is given in SectionlH The first subsection there outlines a general technique for such com- 
putations, while the last two subsections are devoted to the computation of diffusivity for a 
periodic focusing semicircle micro-geometry and related parametric families. 



3.1 Outline of proof of the central limit theorem 

In this section we explain the skeleton of the proof of Theorem |2] as a sequence of lemmas, 
leaving the proofs of the lemmas for Section [5?2] Now Z represents more generally an in- 
tegrable scalar (rather than vector) random variable so that E[|Z|] < oo, having mean and 
slowly varying at infinity; that is, E [(Z^)^] - O(lna). Although Z is more general than before, 
we find it convenient to continue to refer to Z as the displacement. The typical random vari- 
able we wish to apply the below theorems to are the projections of the displacement vector 
Z" considered before. For a channel in bounded by parallel lines, we are mainly interested 
in Z{v) - rvil V2, where v - (vi, V2), with vector (0, 1) being perpendicular to the boundary 
lines. 

Lemma 1. Letij, he the random inter-collision time. Then {ah[a) t)~^Na,t converges tol/Ey[T},] 
almost surely for each t> as a ^ 00. In particular, let 

ah{a) t 



where [x] denotes the integer part of x. Then Na,tl na,t ^ 1 almost surely. 
The proof of Lemma[T]has already been given in the proof of Theorem[T]in Section lZ!2] 

Lemma 2. Let naj be defined as in Lemma\l\ IfY."2'o ^ converges in distribution to 

3^(0, tT)) then so does'E.^^o'^ a'^Zj. 

The next lemma shows that one can work with the a-truncated random variables Z^j in- 
stead of the Zj . 

Lemma 3. The quantityY.'j'^Q ^ [Zj - Zaj) converges to in probability as a— ^00. 

To address the issue of statistical dependence among the displacements, we employ Bern- 
stein's big-small block technique. That is, we break the sum of truncated displacements into 
alternating big and small blocks in such a way that the small blocks are negligible and the 

big blocks are in a sense independent. Let a - 0.01, /3 - 0.6 and define ba,t - ['^f f] and 

Sa,t - [^at]- These are the lengths of the big blocks and small blocks, respectively. Define 
the big blocks as 

Ua,i = ^ a ^a,(i-l){ba,t+Sa,t]+j' 
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for l<i<ka,t, where ka,t is the largest integer / for which [i - 1) [ba,t + Sa,t) + ba,t < riaj - 1. 
Note that ka,t ~ n~_ 
big blocks. That is, 



Note that ka,t ~ n\ f ■ Next define the small blocks as the sums that remain between the 



Va,i - ^ a ^Za,{i-l)[ba.t+s„,t)+ba.t+j' 
7 = 1 

for 1 < J < ka,t, and Va,k„_, = a'^ [Za,(ka,t-mb„,t+s„.t)+ba.t+i + ■■■ + Za,n„,,-l] , SO that 

"o.t-l k„j 
7=0 i=l 

The next lemma shows that it will suffice to consider only the big blocks as the sum of the 
small blocks is negligible in probability. 

Lemma 4. The SMm xfl'i converges to in probability as a — ► oo. 

The proof of Theorem|2]has by now been reduced to showing that zf°i Ua,i ^(0, tD) as 
a—'oo. Therefore, the theorem will be proved if we show that the characteristic function of 
this sum converges to the characteristic function of the normal random variable. By the next 
lemma, the big blocks are asymptotically independent in the sense that the characteristic 
function of their sum can be estimated by the product of their characteristic functions. 

Lemma 5. Convergence |e |exp t^a.ij] ~ nf"i' E [exp (j/i | holds for all fi in R 

as a—- oo. 

Combining the above with the following lemma then gives the proof of Theorem|2l 
Lemma 6. T/ie convergence nf"'jE[exp(/s?7a,,)] ^exp|-yrl>j holds for all s eU as a ^ oo. 

3.2 Proof of the weak invariance principle 

We give now a proof of Theorem [3] and of Proposition HI To show weak convergence in the 
space C[0,oo) := C([0,oo),IR*^) of continuous paths in it suffices to show that the finite di- 
mensional distributions of the projection of X„ along each coordinate axis of U'^ converge 
weakly to those of Brownian motion B and the collection Xa in C[0,oo) is tight (see, for ex- 
ample, (5]). In fact, with regard to the second condition it suffices to show tightness of the 
collection restricted to C[0, f] for all f > (see, for example, |28|). Thus the following two 
propositions are sufficient to prove the theorem. We assume without loss of generality that 
dimension A; = 1. 

Proposition 5. Under the conditions ofTheorem\^ the random vectors {Xa,ti,---,Xa,ti) con- 
verge weakly to iBt^,...,Bt,) as a ^ oo for all I and all ti < ■ ■ ■ < tj e [0, oo) . 

Proof. The proof is by induction on k. For the case fc = 1 we begin by writing 

-1 



Xa,i 



\ 7=0 



+ Ra,i 
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where Raji is the signed distance traveled in the time between collision N„ tj and fi. We claim 
that Ra,ti ^ in probability as a ^ oo. Indeed note that \ Ra,ti I £ \ a~^Zi\f^ | so, for any e > 



P(|i?a,tJ>e)<p(|fl-iZjv„,,J>c)< 



IE[|Z|] 



which goes to as a — ► oo. It follows that Xa,ti converges in distribution to Bt^ by Theorem|2l 
Next we consider the case / > 1. It suffices to show that ^i^a.ti converges in distribu- 
tion to (iBtj as a^oo for any e K'. We first write 

I 

X ^i^a.ti - ilXa,ti + l"^i-2^a,f,_2 + (f'-l + ^/)^a,f/-i + f / (^fl,t/ ~^a,f;_i) 

i=l 

Arguing as in the previous case and using the techniques of truncation and Bernstein's method 
as in the proof of the central limit theorem we conclude that the first / - 1 summands above 
are asymptotically independent from the last so that, by the induction hypothesis, the sum 

ilXa.ti + 1" ^/-2^a,f;_2 + (fi-l + ^fc)^fl,t;-i +il{Xa,t, -^a,t/_i) 

converges in distribution as a ^ oo to 

/ 

1=1 



This concludes the proof of the proposition. 



□ 



Proposition 6. Let t > and define u{a,S) :- sup{|Xa_[, - Xa,v \ :\u- v\<8 and u,ve [0, f]} . 
T7ze?zlimg^olini«— ooP(m(<2.5) > e) = foralle > 0. 

Proof. Let e > and S < t. For simplicity we assume 6 divides t and let n- tIS. The argument 
holds in general with only minor modification. Let = fo < • • • < tn- f be the equidistant 
partition of [0, f] . Observe that 



^[u{a,8) >£■)<! 



max sup Xa,s-Xfl,f 



>e/3 



n-l 



sup \Xa,s-Xa,tj\>el'i 
\s£|tj,f,+i] 



Introducing the notation d[. := Z'-^j- Za,j - ^ Za,t, then for any < j < n - I, the event 

sup||Xa,5-Xa,f^. I :sE [fj,fy+i]| > c/3 implies that maxN^^^^,<:jc<:N„,,j^^+i \Dq\ > ael6. Next let 
ns = [ "l^jfjf ] and let A{j,S) denote the event Na,tj^-^ - Na,tj + I < «5. Then the probability 
p(sup^£|j^._(^.^j |^a,s-^a,iy| >e/3,A0',5)j is bounded above by 



max 

^Ar„,f.<fc<jv„,,^.^^+ 



\d^\ > ae!6,AU,S) 



<2P 



max 



k-1 

Lzi 

i=0 



\Dq\ > ae/6 



max 

l<fc<)!5 



> ae/6 



<4P 



ng-l 
i=0 



> ae/6 , 
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where the last two inequalities follow as in the proof of Lemma[2l Therefore, 



\S£ltj,tj + l] 





n-l 


1 




n«-l 






>c/3 










> ae/6 












i=0 









1 


1 


ns-l 


\ 




t 


/■ 




ng-X 






n 








> ae/6 




4P 






> ael6 










i=0 












i=0 







Notice that the number of collisions in the interval \tj, is precisely Na,tj+i - Na,tj + 1 and 
so by Lemma[T] {Na,tj+i - Na,tj + ij ^ 1 almost surely as a — ► oo. From this it follows that 
P[A{j,d)'^) —- almost surely as a — ► oo, independent of j and 5. From Theorem|2] 



1 



ng-l 
i=0 



>el6 



as a ^ c». Letting (5 ^ on the right-hand side above then gives the result. 



□ 



Proof of Proposition^^ The equality in the statement of theproposition is a standard fact on 
the mean exit time of Brownian motion from an interval. We prove here only the convergence 
of mean exit times. Note that by the continuous mapping theorem T{Xa) converges in distri- 
bution to T {B) . (Of course t is not continuous on all of C [0, oo) but it's not difficult to show that 
it is Pg -a.s. continuous.) Therefore, to show the convergence of mean exit times it suffices to 
show that the collection of T{Xa] is uniformly integrable. That is, it suffices to show that for 
any e > there exists M > such that EfJ [t1{t>m) ] < e for all a. 

Note that for all e > there exists 8 e (0, 1) such that Pf (t < e) > 5 for all xe{-L,L). Since 
T [Xa) converges to t {B) in distribution it follows, similarly, that for any c> there exists some 
5 £ (0, 1) and ao such that PJ(t < e) > 5 for all a> oq and for all x £ {-L,L). If we let e = 1 and 
let 5 be the corresponding value in (0, 1) then it follows by induction and the strong Markov 
property that Pq (t > fc) < (1 - 5]'^ for every positive integer k and a > oq. Therefore, if we 
choose M' large enough so that 'Z.'l?i^,{k+ 1)(1 - 5)*^ < e, then for a> oq 



>M'] 



k=M' 



< e. 



k=M- 



It is also straightforward to see that there exists M" > such that EfJ [Tlj^>jvf")] < ^ for <^ < "^o- 
Letting M = max{M', M"] then gives the uniform integrability. □ 



4 Examples 

This section is devoted to showing how the diffusivity D encodes surface microscopic struc- 
ture when our operator P represents a random reflection. The structure of the section is as fol- 
lows. The first subsection gives a general outline for computing D independent of any given 
surface microscopic structure. The second subsection computes D in the case that the sur- 
face of the walls is given by a periodic arrangement of focusing semicircles. The last subsec- 
tion gives D for certain parametric families of surfaces derived from the semicircle example 
of the previous subsection. 
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4.1 General Technique 

While in general, under Assumption[T] the diffusivity is given by 



fD = lim E 

a— ♦cxD 



7=0 



Ha] 



where I(,a) is the interval given in Theorem|2l it is possible to consider a significantly reduced 
truncation without altering the value of D. 

Let T] e (0, 1) and define Jia) :- {x : exp (log'' a) < \x\ < al log^ a}. The following proposition 
shows that we may use the truncated displacements Z^'''' in computing D so that in fact a 
vanishingly small cone of trajectories determine the diffusivity. 

Proposition 7. Under the assumptions of Theorem\^ 



tV = lim E 

a— 'OO 



7=0 



/(a) 



Proof. As in the proof of Theorem[l] the second moment on the right above may be expressed 
as the spectral integral 



N„,t-l 
7=0 



Ka] 



1 l + AE[Na.t + Oa]] (h{a) 



ah{a) 



\ a 



Z 



Ha)\ 



Uz]l.a){dX). 



where H^/m [dX] = ||Z-»'" || (Z^''",n(<iA)Z'f''"> is a probability measure on the spectrum of 
P. Now observe that if W Z^''"^'^'''^ then 

(z-^f'^ntdAjz^'"') = (z^"'' - w,n[dX) [z'^"^ - w]) 

= (Z'^"^ , mdrtZ"-"^ )-2(W,n {dDZ"-"^ ) + {W,mdA]W) 
and \\W\\l^ 0{\n^ a). Therefore, 



Z 



Jia) II 2 
2 



^7(fl)||2| 
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;2||PV||2||Z^'«'|i2 



+ ||W||? = 0an'' a) 



for some rj' e (0, 1) . Multiplying both sides of the inequality hy h{a)/ a - 1 / In a and taking the 
limit as a — ► oo implies that Z^^"^ and Z^^"^ grow at the same rate, and that in the limit formula 
for D we can use Z^'"^ rather than Z^'"'. □ 

With the previous proposition in mind we use the shorthand notation Za ■- Z-'^*"' through 
the rest of the subsection without the risk of ambiguity. 
If we expand 



JV„. 



-1 



a ^ ^ Za j 

7=0 



\2 




Na.,-1 








= E 


E 4, 

7=0 


+ E 


l<i<j<]Va,,-l 



the limit of the first term has been shown in the proof of Theorem [T] to be Dq, which is inde- 
pendent of the microstructure. Let C{a) be a function that increases slower than ah{a) but is 
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otherwise to be determined by a specific microstructure. For tlie second term we break the 
sum into two pieces as follows 



2 L 

l<i<j<N„j-\ 



Q<\i-j\<C{a) 



E 

C(a)<|i-j|<W„,,-l 



The first term on the right above wiU be determined by the microstructure, but the second 
term actually vanishes. 

Lemma 7. Wehave\ima^ooa~^^['Lcia)<\i-j\<Na,t-i^a,iZa,i] =0. 
Proof. Observe that 

Nr,,t-C(a) N„,,-C{a) 

L \E[ZajZa,j]\^a-^ ^ ^ |E [Z„,iZ«j+C(a)] | 

C(fl)<|i-j|<Ar„,,-l i=l ]=i 

Na,,-C(a)N„,,-C(a) 

= 2fl-2E[z2o] X E |corr(Z„,,-,Z«j+c(a))| 

iV„,,-C(a)W„,f-C(a) 

<2a-^E[zl^] L \piC(a] + j-i]\ 

i=l j=i 
A/f N„,,-C(a)N„,t-C{a) 

!=1 ;=i 



ah{a) 



where M > is a constant and < p < 1 is the essential spectral radius. Further, 



Na.t-C{a)N„j-C{a) C(a) 
^C{a)+j-i _ P 

1-p 



E E 



l_ „Na,t-C(a)+l 
Na,t - C{a] + 1 



1-p 



Taking expectation and letting a — ► oo then gives the result. 



□ 



4.2 Semicircle Microstructure 

This section is devoted to computing the diffusivity of the channel system whose walls consist 
of a periodic focusing semicircle microscopic structure. The first subsection is devoted to a 
thorough analysis of the operator P for this geometry. We show that P is quasicompact and 
give an explicit formula for P for a certain range of pre-coUision angles. The second subsec- 
tion is devoted to the computation of CD using the method outlined in the above section. 



4.2.1 A CLOSER LOOK AT _P 

To show that P is quasicompact we employ the technique of conditioning (see 1161 and [17] for 
more details). The general idea of conditioning in this setting is to obtain a compact operator 
by considering P conditional on the event that trajectories satisfy a given property. As we will 
see, the compactness of this conditional operator will imply P is quasicompact. 
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Let Q denote the the microscopic cell bounded by the semicircle and its diameter such that 
= To U Fi where we renormalize so that the diameter Fq is identified with [0, 1] and Fi is the 
semicircle with radius 1/2. Let ^gtr) e [0,:^] denote the angle between the outgoing vector V 
and Fo given that the trajectory enters Q with angle 6 at position r. 

We define a measurable partition of M = / x y = [0, 1] x [0,n] as follows: let Mi c M be the 
subset of initial conditions whose billiard trajectories undergo exactly one or two collisions 
with dQ before returning to Fq. Define Pi as P conditional on the event M\, and similarly 
define P2 from M2-M\M\. More precisely, if we let Mj[6) - [r e I : [r,6) e Mj} and define 
aj [9) = A[Mj [9]) for each 7 = 1, 2 and e then for each / e L°"[V, n), define 



[PjfW) 



0, aj(9]^0. 



We call Pj the conditional operators associated to the partition Mj . Note that it makes sense to 
write Pf - aiPif + 1x2 ^*2/- Let /ij be the measure on V defined such that d^j - ■fXT^i^jg^'^A'- 
It follows that Pj is self adjoint on I? { V, j^j ) . 

Next note that Pi has an integral kernel. Let - {r e I : {r,6] e Mi]. Then Wg is the 
countable (or finite) union of open intervals Wo j for which the restriction ^'gj - ^g\wgj is 
a diffeomorphism from Wgj onto its image Vgj. Define Tgiq)) - Y.jXVgj{^)f^B,iW)~^: where 
Ae,j{(p) = \ sin(/). Let 

^ (Ax^)(Mi)Fe((p) 

m{9,q)) = — — — — . 

ai[9]ai{q)) 

It follows by way of change of variables that the operator Pi on L^{V, fii) is given by 

(Pi/)(0)= f f[ip)a)i{9,(p)dni[(p). 
Jv 

To show Pi is compact, it will suffice to show wi is square integrable onV xV. To this end, we 
first look at the function I'g (r). By the symmetry of the semicircle I'g (/ ) satisfies 

Wg{r)^7i-W„.ga-r). 

It follows that it suffices to consider only 9 e {0,7i/2). 

Proposition 8. The function ^'g (r) has the following properties: 

1. Let9e {0,ji/2) andletn> 1 be the number of collisions a trajectory with initial data [r, 9) 
makes with the semicircle. Then 



2nsm-^[{2r-V,sm9) + mi-9, re [0,1/2] 

2nsm-^{[2r-\)sm9)-[n-2)n-9, r£ (1/2,1]. 



2. Let9 e {n/4,n/2). Then ^g{r) has the following points of discontinuity in [0, 1]; 

(„)_1 sm[{mi-9)l{2n+\)\ („) _ 1 sin[((77- l);r + 0)/(2«+ 1)] 
^° ~2 2 sine ' ~ 2 2sin0 ' 

n>\. 
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3. Let e [0,71/4:). Then (r) has only one point of discontinuity given by 

, 1 sin0/3 

r -- + . 

2 2sin0 

We remark that if e {Q,nlA) then at most two collisions in the semicircle are possible and 
hence the formula in 1 . is only valid for n- 1,2. Further, we make more precise what is meant 
by the points of discontinuity given in 2. and 3. If r e < '"{^'j' then the initial conditions 
(r, 9) give a billiard trajectory which makes only one intermediary collision. And for « > 2, if 
r e ^rg"', ^g""^'! U ['"{"~^\ r{"'j, then the initial conditions (r,0) give a billiard trajectory which 
makes n intermediary collisions. For Q < 6 < nIA, the situation is simpler. The initial con- 
ditions (r,0) for r e (0, r'), give one intermediary collision, and for r in the complementary 
subinterval of I, there are two intermediary collisions. The proof of the proposition is by ele- 
mentary trigonometry. 




Figure 4.1: Definition of r. 



Propositions. The operator Pi onI?{V,pi) is compact. 

Proof. As Pi is given as an integral operator {P\f){9] - Jy f{<p)CL>\{6,ip)dp\{q)) it suffices to 
show that wi e L^iV ^V,ii\^ p\). By symmetry, we have the identity - (i)\{n-6,n- ip) 

and the square of the L^-norm of wi is 

nnl2 r 

2 i i (D\[B,q))dpi[q))dpi{e). 
Jo Jv 

Thus we show that the following two integrals are finite: 



C4.1) / / (u\{ff,(p)dp.i[(p)dpi{e), \ i (n\(B,q))dpi[(p)dni 

Jo Jv Jnii Jv 



m. 



First consider the second integral. Let 7i/4 < 9 < n/2. Using the notation established in the 
previous proposition and further above, Wg is given by the single interval Wg i - (rg^',rj^'). 
Correspondingly Vq^i = {7il3-BI3,7i-B/3]. Moreover, ^e,i(r) = :T:-0 + 2sin"i((2r- l)sin0). 
It follows that 

, , 4sin0 
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Now observe that 

J71I4 Jv Jjiii Jv ai[e)ai[(p) 

rnl2 



~ 1/4 Jvg^, ai(0)o:i(^)|^;j(^-i(^))|2sin(^ 



Now because 6, (p are bounded away from zero in the above integrals and because ai (0) > 
for all 6, it follows that the above integral is finite. 

We show that the first integral in 1 14.1) is finite. Let < 6 < nl4. Here Wg is given by the single 
intervalWg,! = (0, r') and Vg,! = [71-36,71-913). Moreover, ^e,i(r) = :?r-0+2sin-i((2r-l)sin0) 
as in the previous case. Hence \^'gCi'g\((p))\ is also as above. It follows that it suffices to show 
the following integral is finite: 

r7ii4 /-jr-e/a sin^ l^^l 

^ '-d(pde. 



Jo Jn- 



-36 sinfc/sin(^ 

Note that for 71 - 39 < qx 71 - 9 1 3, we have (p>6, which implies sin j £ sin^ (p. Therefore, 

Jrnii /-Tt-e/a sin^ l^i^l pnl4 rn-813 a[ri(n 

/ 5^ '-dipd9< \ / — -d(pd9 
Jji-se sin sin Jo Jn-w sin0 

C siuM 

= / / 30 -dwdO 

Jo J0I3 sin0 

r"''^ 1 f (30)2 -(0/3)2 

~ Jo sin0 



d9, 



which is finite. □ 

Having shown that Pi is compact we are now ready to show that P is quasicompact. This 
this a consequence of the following general fact. (See Theorem 9.9 in [27 1 .) 

Proposition 10. LetK and T be bounded self adjoint operators on a Hilbert space and suppose 
thatK is compact. Then the essential spectrum ofT + K is contained in the essential spectrum 
ofT. In particular, if\\T+ K\\ - 1 and \\T\\ < 1, then the spectral gap jiT + K] ofT + K satisfies 

j{T+K)> min{l - || r|| ,r{K)}. 

The quasicompactness of P then follows from letting K - aiPi and T - azPi and noting 
thatl-liril >infai>0. 

4.2.2 Computation OF D 

In this subsection we use the shorthand Za :- Z^^"^ and Zaj :- Z^'^^Oj) where J{a) is as 
given at the start of Section fOI 
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We show, for the example of Figure 12.21 that D - ^jz^ where ( - -|log3. Following the 
discussion in Section fOl we aim to compute 

lim X E[Za,iZa,j] 

"^'^ 0<\i-j\<C{a] 

where we choose C{a) - \og^ log a. By stationarity, we are interested in computing E [Za,oZa,j] 
for j < C{a). That is, we are interested in the first C{a) collisions of trajectories with a shal- 
low initial angle. Now although we've chosen these truncations of the between collision dis- 
placements because, as seen in the previous section, the transition probability kernel has a 
straightforward explicit formula for such pre-coUision angles, they are not without their own 
complications. It's clear that to keep track of the trajectories whose 7th displacement falls out 
of the truncation range, making Za,oZa,j vanish, quickly becomes intractable. For this rea- 
son, we introduce the following so-called widened truncation, which, as we will show, will not 
change the diffusivity. 

Consider a trajectory for which Zq e J{a). I follows that, for large enough a, 

\z(?.j&Q)\<\Zj\<\z(?.-i©Q]\. 

for all i < C{a). A straightforward estimate then shows that 

3"-'' expdog'' a)-Ci<\Zj\< 3^ al log^ a+Ci 

where Ci > 0. 

Define K{a) := [x : S"^'''' expQog'' a) < \x\ < 3^''''a/logT' a + Q}. For the rest of the subsec- 
tion we will consider the new truncated displacement z^'"' as well as Z^*''', which we will 

continue to denote Za- Note that each zf - Z^''''(0;) is nonzero. Moreover, the following 

J ^ 

lemma shows that it suffices to compute E ^ZajZ^''"^^. 
Lemma 8. The following equality of limits holds: 



lim 

a— 'CO 



X a-\[Za,iZa,j]^ lim a-^ ^ E[z„,iZp'] 



0<\i-j\<C{a) 



0<\i-i\<C(a) 



Proof. Define h := {x : 3 '^'''^ expQog'' a) < \x\ < exp(log''fl)} and I2 :- [x : a/log^a < |x| < 
3'^^"^ a/logf a+ Ci). Note that K{a) - Ii[a) u hia)- Hence it suffices to show 



lim a 



- E ElZ, 

0<\i-j\<C[a) 



,,;(ZJ1'''' + Zj^'^')] =0. 



Observe that 



E 



< a-V 



-111 



/i(a)V 



-HE 



1/2 



O (a"^ log"^ adoglog a) ^'2) , 



where the last step is due to E | Z^ ^ j - 0(log a) and E \zb '''' j 
contains roughly 0{Na,tC{a)) - O^ " such terms, the result follows. 



- O(logloga). Since the sum 



□ 
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Let qi - Zj/Zj-i for i>l. The next lemma is a key technical tool in the computation. 

Lemma9. Letl< 7<O(|log0|). ThenE[qi--- qj \Qo^ 6] ^ (;i + ABi'^B'^ + 0[B'^) for constants 
A,B independent of 6 . 

Proof. Observe that for 6 sufficiently small, using the integral kernel for P derived in the pre- 
vious subsection 



^ / nii-di3 COS ( '^"^2 " ] r^s cosj^^J \ 

E\qi\@o-9]- / cotw -dw+\ cotw dw\ 

" ^ CQtByL-3B ^ 4sin0 ^ je/3 ^ 8sin0 

1 f 1 {<p-e\ ^ 1 (<p+e\ ^ ] 

— / cotMcos \dw+—l cotcpcos \d(p\. 

cosei AJeis ^ [ 2 j ^ 8je/3 ^ I 4 j 

One may check that 

[ cotmcosf^^] <i«) = 21og3 + Di0^ + O(0*), 
cosBJeis \ 2 I 

where Di is a constant, and likewise for the second integral above, albeit with a constant 
different from Di . The case j - I then follows. 
Let Ve = (0/3,30) u (:?r- 30, TT- 0/3). Observe that 



£[^71^72 I ©0 = 0] =E[E[^7i^?2 I ©0 = 0,01 = <p] 100 = 0] = [ E[qiq2\&o^d,&i^<p]PiO,dcp) 

JVg 

= [ ^E[q2\&o^e,&i^(p]P[6,d(p)^ [ ^[( + A(p^ + Oi(p^)]P{B,dq)] 

JVg Z{B) JVg Z[U) 

= CE[t7il0o = 0]+ [ 4^(V + O(/))P(0,rf(p) 

JVg Z{B) 

^e + Ae^ + 0[B']+ f ^[Aip' + 0[cp'))PiB,dcp). 

JVg Z(0) 

Next, one may check that 

[ ^[A(p^ + 0{cp'^))PiB,dcp)^-^i-l [ cot(p{A(p^ + 0[ip^))cosi^^]d(p 
Jvg z{0) ^ COS0 I 4J0/3 [ 2 j 

+ - f cot(p{A(p^ + 0{(p^)]cosi^(^^]d(p\^ABB^ + 0{B^) 
8Je/3 14^ I 

The case j - 2 follows. The rest of the argument follows by a similar induction argument. □ 

With the above lemma in place we are ready to compute the correlations. 

Lemma 10. For 1 < j < C[a), E [z„,oZj*^'°'] = 4r^(J A{a) + ABj~^r{a), where Ma) ~ loga, 
T{a] - 0(exp(-21og'' a)), andA,B are the constants given in Lemma\^ 
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Proof. Recall the interval J{a) - {x : exp (log'' a) < \x\ < a/log^ a}. Observe that 

= 2 f Z{efE[qi---qj\&o^e]l\na)^iid6) 
Jo 

= 2 ZiOf [(j + AB^-^e^ + 0(0*]) ^{dB). 

We remark that because v^^e are only considering here such that exp (log'' d) < Z{6)< a/log^ a, 
and hence | log0| < Clog a for some constant C, it foUov^rs that j < C{a) is sufficiently small so 
that we may apply Lemma|9l 

It is straightforv^^ard to compute Jq Z{d)^'ij{a] ii{d6) ~2r^ log a. And moreover 

f Z(0)2o(02)l/(a) tiidB) = r O(0)l/(„) dB^O ((expdog'' aW^] . 
Jo Jo 



The result now follows. 

The summation of correlations is the final piece to our computation. 
Proposition 11. LetU{a) ^j:o<u-j\<c(a)^{Za.iZf'-'''>j. Then\ima^ooa-^U{a)= 
Proof. Observe that 

C(a)-1 r 1 C(a)-2C(a)-!-l 

(Na,t-C{a) + 1) J [4r^(jA[a) + ABj-^na)) 

C(fl)-2C(a)-i-l 

L E 

i=i j=i 



□ 



U[a) = 2 



= 2 



C(a)-1 C(a)-2C(a)-!-l 



! = 1 

C(fl)-1 

E 



C(fl)-2L(a)-i-l 

+ E E (4rVA(a) + ^5-'"^r(fl) 



= 2 



.rr nr . , B«fl) + ^^^M - 2) - ^2 (CM - 1) > 

^r(fl) (Ar„ f - C(a] + 1) .. + ■ 



B[B-l) 

rCla) . 



B(B-1)2 



2 f C-^"! +([C[a)-2)-(''[C[a)-l)' 

+ 4r^A[a)\[Na.t-C[a) + i y , + 



C-1 



Dividing by a and letting a ^ oo gives the result. 



□ 



The value of D then follows from adding the value from the proposition above to 4trv In as 
discussed in Section l4?T] 



4.3 Parametric Families 

In this section we consider three different parametric families which are derived from the 
semicircle. Of primary interest will be how the diffusivity of the limiting process for each 
family changes as a function of the parameter. 
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We begin with the family formed by adding a middle wall of height h to the semicircle as 
shown in the figure above. Suppose h < 1/2; that is, the wall does not extend to the center 
of the semicircle. It is apparent by inspection that trajectories with a sufficiently small pre- 
collision angle will never, so to speak, notice the middle wall. And moreover, those trajectories 
with initial data (r,0) that do notice the wall will behave like trajectories in the semicircle 
with no wall with initial data (1 - - 6) by symmetry. That is, if ^g{r) denotes the post- 
collision angle of a trajectory with initial data (r, 6) in the middle wall geometry and similarly 
= *I'g(r),then^g(r) - ^;j_e(l-r) for initial conditions (7",0) for which the trajectory hits 
the middle wall. It follows that the operator P;, for the middle wall geometry is quasicompact. 
Moreover, as pointed out earlier, the diffusivity D;, for P;, depends only on trajectories with 
arbitrarily shallow angles by the formula in Proposition [T] It follows that the diffusivity D{h) 
is constant and equal to the diffusivity for the semicircle with no wall for all h< 1/2. 

Further, using the symmetry 'l'g(r) - ^j,-0(l-r] and following the proof given for the semi- 
circle, the diffusivity D (1/2) for the geometry with a middle that extends exactly to the center 
of the semicircle is given by 11(1/2) = ^ where ( = -^logS as in the semicircle. We sum- 
marize these facts as follows. 

Proposition 12. For the middle wall modification of the semicircle with middle wall height h, 
Dih) - D, for h< 1/2, where D is the diffusivity for the semicircle with no wall, and for h- 1/2 
the diffiisivity is given by 

D(l/2) = -4, 

where(= -l/41og3. 

It is interesting to note that the diffusivity is not a continuous function of the parameter. 
We also remark that the case h > 1/2 when the middle wall extends outside the semicircle 
requires a different analysis altogether, which we leave for a future paper. 

Next we look to the geometry formed by splitting the semicircle and adding a flat bottom of 
length h e (0, 1) as shown above. We renormalize the size of the semicircles so that they have 
radius (1 - h]/2. We also establish the notation a = {1 + h)/2, b- {1- h)/2. While qualitatively 
similar to the semicircle, the angle function ¥g (r) requires a new detailed analysis, which we 
sketch here. 

We begin by noting that by symmetry it again suffices to consider only pre-coUision angles 
9 e [Q,n/2). And moreover, for sufficiently small pre-coUision angles 6 at most two interme- 
diary collisions are possible within the cell. In the discussion that follows we consider only 
such 6. Let r' e (0, 1) be the point of entry for which ^g[r) is discontinuous. That is, r' is 
chosen such that trajectories with initial data (r,0) for r e (0, r') experience one intermedi- 
ary collision, and for those with r £ (r', 1), there are two intermediary collisions. It follows 
that T'^(r) = 7i-B-2pi(r) for r e (0,r') where Pi = j3i(r) satisfies fosin;6i = {a-r)smB. We 
may also characterize r' as the value of r that satisfies asm{6 - 2p\) - fosinjSi. From these 
observations it follows that 

3 + h 1 — h 

"¥^0)^71 6 + 0(0^], lim Whr)^ji 0+O(d^]. 

" l-h ^ ' r^ir']- " 3+h ^ ' 



Following the notation established in the discussion on the semicircle. 



3+h , c,, l-h , 
e + 0{B'),n--—9 + 0{B') 



l-h ' ' 3+h 
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Moreover, if we let B = I'g (r) it follows from implicit differentiation that 

1 4 sine 



&'{r) ■■ 



l-h, 



,1 0+0-7l \ 



which is the corresponding value in the semicircle case except for the factor of{l-h) ^ . 




Figure 4.2: 



For r £ (r',1), we have "Vgir) = 2jSi (r) + 2/32(r) - 0, where Pi{r) is as above and P2 = Pzir) 
satisfies bsmfiz - ^sin(2/3i - 9] + smjSi. It follows by symmetry that 



3 + h 



e + o[B^],^—^B + o[e^] 



Further, by implicit differentiation again 

©'{r): 



l + h 8sin0 
a-h)^ cos(^)' 



which again resembles the semicircle case but for the factor of [l + h){l - h) ^. As in the 
computation of the diffusivity for the semicircle we find the following. 

Proposition 13. For the flat bottom of length h e (0, 1) modification of the semicircle, we have 
the following value of diffusivity as a function ofh: 



D{h)^ 



4rv l + (k 
Tt i-(h' 



where 



l + 3hl-h 3+h 
■log- 



4 l + h 1-h 

Notice that the limiting case h - gives the diffusivity of the semicircle while h - 1 gives 
(h-0 which gives the diffusivity in the case that at each post-collision angle is chosen inde- 
pendently according to the distribution /i. 

The final family of interest is formed by adding a flat side of length h e (0, 1) between semi- 
circles as shown above. The key observation here is that the operator Pf^ corresponding to 
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such a geometry can be thought of as a sum of conditional operators given by Pf^ - {l-h)P+ hi 
where P is the operator corresponding to the semicircle geometry and I is the identity oper- 
ator. Such an operator is a generalized Maxwell-Smoluchowski model as discussed in the 
introduction. Next notice that since P and I commute the spectra of P/j and P are the same, 
and in fact every spectral value of P;, is given by (1 - /i)A + ^ for some unique spectral value A 
of P. It foUoviTs that 

where is the spectral measure as given in the statement of Theorem|2] and similary Yl'j, is 
derived from the projection valued measure n'' associated to Pj, given by the spectral the- 
orem. The following proposition then follows from the discussion above and the spectral 
formulation of the diffusivity. 

Proposition 14. For the flat top of length h e (0, 1) modification of the semicircle, we have the 
following value of diffusivity as a function ofh: 

^ 1 ^ h 4rv 4rv ( l + ( 
T>{h]^ -D + ; = + 



l-h \-h n n{l-h]\ 1-0 
where 2) is the diffusivity for the unmodified semicircle and( = -1/4 log 3. 

5 Additional PROOFS 

We collect here proofs of some of the more technical propositions and lemmas from earlier 
parts of the paper. 

5.1 Sketch of proof of Proposition [2] 

The proposition results from tedious but elementary and straightforward integrations. We 
show a few steps to convey the flavor. 

By identifying the tangent space to G at q e dG with ffi the unit normal vector n is 
identified with e„ - (0, . . . , 0, 1) and is identified with the half-space consisting of vectors 
V = [vi, V2) e ffi IR""*^ such that {v,en) > 0. For such a post-collision velocity v, the point 
of next collision with the channel boundary is q + TV, where ij, -2r{v, e„>/| ^2!^ is the time 
interval between the two collisions and Z{v) -2r(v, Cn) vi/l ^21^ is the between-collisions dis- 
placement vector in "horizontal" factor K*^ as represented in Figure lzTl We assume that v has 
the following general form: 

dv{v) = C{v, n)f[\ vf-] dV{v], 

in which /(x) is a nonnegative function on [0,oo) such that x" f [x^] dx < 00. We allow / 
to be distributional (a delta-measure concentrated on a fixed value of | v\) so as to include the 
case where v is supported on an hemisphere. Let S+ :- S"~^ be the unit hemisphere in H. 
Then IHl is diffeomorphic to S+ x (0,oo) under polar coordinates v ^ {w,p), where w - v/\ v\ 
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and p - \ v\, and dV{v) - p" ^ dV^^t^iw) dp. Also define the notations S+a '■- S+ n ll-ll(fl) and 
E[a) -.^ Ey\\Z\^ ,U{a)\.Then 



/ p"f[p'^)\Z{w)f-{w,n)dV,^^[w)dp = C'\ \Z{w)\^{w,n) dV,^, 

'^S+a •'S+a 



iu>), 



where we have used that Z(,v) - Z{w). Here and below, C, C' , C" are positive constants that 
can be obtained explicitly. 

The image of S+a under the projection map n : w e S+ ^Uj e U"~^ contains the ball B^~^ 
of radius (l + {2/a)^] , and it is equal to B^"^ when n-k-1. The volume elements rfKph 
on the hemisphere and dV„-\ on the ball are related under nhy {w, n)dV,^i,{w) - dVn-i(w). 
Also let xi - 7ii(x) be the natural projection from the unit ball in to the unit ball B*^ in 
the horizontal factor r'', and X2 - nzix) the projection to the complementary factor IR"~*^~^ in 
R""^. Then the limit as a ^ oo of /^^^ \Z{w)\^{w,n) rfKph(w) is, up to multiplicative positive 
constant, the same as the limit of 

r (l-W-b^^)^ f r Ixil^(l-|xi|2-|x2|2) 

JBr' (l-|Xl|2)^ M%^f-' (l-|Xl|2)2 

2Vol(B 



— — /|Xi|2 l-|xi|2 ^ dVfc(Xi) 

71-fc+l Jb* 

—Vol (b"-*^-!) Vol (b*^) shl-s)'^ds 



where the iterated integrals were carried out in polar coordinates, in which volume elements 
are related by dV,i{v) - \ v\"~^ dV,j,i,{v/\ v\) d\ v\. The limit of the remaining integral, as a goes 
to oo, is a Beta-function of the exponents of s in the integrand; it, and the volumes of unit 
balls, can be written in terms of Gamma-functions and further simplified. 

The expected values of t i,{u) - 2r (v, n) / [vzl^ are shown by similar computations to take the 
form 

2j- r { """"^ l 1 r f— ] 

E^^[Tb] = —-^^^nn,k], E^^ [t^] = r(2)SM)^ -l^l[n,k), 

where s is the speed (or radius) of the hemisphere on which p is supported and 

{x,ri) Tf^ 
I[n,k):^ ——dV^,^{x)^ 



{n,k):^ [ 

JS-^ 



|X2|2 (n-fc)r(f)' 



S+ being the hemisphere of radius 1 in H. 

Finally, observed that Ey |(Z^)^ j - Ey [\Za\^] Ik since v is rotationaUy symmetric in the [ 
subspace. The proposition is now a consequence of these observations. 



5.2 Proofs of CLT lemmas 

-1 ^ ^. ^N„,-\ 



Proof of Lemma]^ WeagainemploythenotationD^ := Za,t. Since Za, 

^"=o~^ ^a.j + Dq"'' it suffices to show D^"' ^ in probability. Before getting to this directly, 
we begin with a few technical remarks to be used in what follows. 
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Let V > 1 and define 



"1 = 



l- — \na,t 



+ 1, «2 = 



l+ — \na,t 



Then ¥{n\ < Na,t £ ?^2) > 1 - e for a large enough. Let a* -e^a, n*^- ria'j. Observe that 



n^ah[a)t a*h{a*)t 
e Ua t £ e < <n„f + \. 

E[T6] E[Tb] 



Therefore, 112 - ni< e^^na,t ~^-n*af ^'^^ now ready to address the convergence in prob- 
ability. Observe that 



Df''\> ae]^f^v[Na,t^ i,\Df-'\> ae] 



"2 



]=ni 



Notice for the second term above that 



E P [Na.t = I > ae) + L P (iVa,t = I > «c) 

jyi«i,n2l 



J?'[ni,n2l ;Y|ni,n2l 



For the first term observe that 

"2 , I , "2 



And notice that 



y = > ae] < y PuVa.f = 7, max |Df|>ae<IP max lof I > <2C 



max Dj > 

\na,t<k<.n2 1 I 



Similarly, 



It then follows that 



max 

nfi,t<k<n2 



max 

l<fc<n2-nfl,i 



max Pj > fle| < [ 



ni<fc<n2 ' 



t 


\ 


y, Za,i 


> ae 


i=na,t+l 





k 


\ 




/ 


k 


\ 


> ae 


< P 


max 




> ae 


i=l 






l<t<«2-ni 


1=1 





max 

l<fc<n2-ni 



k 



> ae 



max pi >ac|<2P 



max 

l<S;<n2-"i 



y. -^u,; 



<2P 



max 



it 

y. -^fl,; 
i=l 



> ae 



Next observe that by the symmetry of Z^,/, 

Za,i\ > ac,maxi<s-<„.^ > aej 

IP ^maxi<j;<„* ^ 



max 

l<fc<;!*, . 



fc 

X/ Za,i 



> ae 



1 

> — 

2' 
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which implies 



max 

isfcsn* 



k 



> ae 



<2P 



> ae, max 

l<fc<n* , 



fc 



> ae 



<2P 



> ae 



We siammarize and observe that by assumption 



722 

E P(Na,t^j,\Df"'\>aej<4P 



j = ni 



E 



1=1 



> ae 



:4P 



E 



8 r"' 

71 tV 



as ^ oo. In letting e ^ the last line above vanishes. □ 

Proof of Lemma\^ Let hia) :- {x : a/log^(a) < \x\ < alogloga} and I2 {x : \x\ > alogloga}. 
It will suffice to show that each of the sums a'^zj'^"^ for i - 1,2 converges to zero in 

probability. 

We first consider the truncation by Iz- As there is no ambiguity, we use the shorthand Za := 
Z'^ and Za,j = Z''^ {©j). Let e > and observe that 



1 








'Na.t-l 




00 


' ^ 1 






E 


> ae 


< P 


E \Za,j 


> ae 


= E p 


^ |Z„j| > ae 






7=0 










]V=1 


U=o ] 





- E 'P'l-^; ^ ^2 for some 7,0 < 7 < l] P{Na,t = ^ 

JV=1 
00 N-\ 

]V=1 j=0 

= E[]V„,j]P(Zoe/2). 

Next note that 

Jrcot"' (alogloga/(2r)) 
sinSde 


= 1 - cos (cot"'^ ((2logloga/(2r))) - 0((aloglogfl)~^) . 

Moreover, a straightforward application of Lemma [1] and the dominated convergence theo- 
rem shows that E [W^.t] - 0{ah{a)). It follows that 

1 





iVfl,f-l 








E Zaj 


> ae 













llogfl- (logloga)^ 



as ^ CXD. 

Next we consider the truncation by I\ . As before we use the shorthand Za Z'^ . To show 
that the sum "q a~^Zaj converges to zero in probability, we use Chebyshev's inequality. 
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To this end, observe that as in the proof of Theorem[T] 



1 l + XE[Na,t + Om] h{a) 



ah{a) 



ZafUaidX), 



where we recall Da - \\Za\\~^ {Za,Yl{dX)Za) is the spectral measure associated to Za- Note 
IZal^ - OQoglogfl) and^(a)/fl= 1/ log a, while all other factors in the integrand are bounded 
as — ► oo. The result follows. □ 

Proof of Lemma\^ Because of the slight difference in definitions, we show separately that 
^ V^a.i ~^ and , ^ 0. To prove each of these, note that by Chebyshev's inequality, it 
suffices to show 



E 



■0' ^KkJ 



•0, 



respectively. We start with the first. Observe that 



and 



\ '=1 



E a'^Zaj 



l<i<;'<fcfl,f-l 



Ja 1 



1 + A 5a,t + 0(l) hja) 
■X ah[a) a 



\Zafna{dX)^0{iahia)f-^) 



Therefore, fc^^E | V^j j - O [[ah{a))" ^] 0. The number of terms in Zi<i<;<r„ t-i Va.iVaj is 

(fcfl.t - l)(fca,t - 2) and each term VajVaj itself contains ^ terms of the form a~^Za,mZa,n- It 
follows that 

2 E E[Va,iVa,j]^O{[ah(a)f^'-P^^^''-')^0. 

l<i<7<fco,,-l 

Finally, the sum in , contains by definition less than ba,t + Sa,t terms. Just as above, it 
follows that Ef 1^2 ] < Oi[ah[a))P-'^) ^ 0. □ 

Proof of Lemma\^ Let Wi = expiifiUaj) for 1 < / < ka,t- Also let a^^,. - E [ w^] - E^ [w ,]. Define 
£(u;i, if2) := E[u;i wz] - E[wi]ElW2] . Because of the small block gap of size Sa,t between big 
blocks, |£(wi. If 2) I = Ictu^i ||(Tu;2l|corr(u;i, W2] \ < 4|corr(wi, W2) \ < 4\p{Sa,t]\- Moreover, 

\E[WiW2W3] -E[Wi]E[W2]ElW3] \<\E[Wi,W2W3)\ + \E[Wi]E[W2,W3)\ 

= \a,i,i\\(Tw2W3\\C0n(.Wi, W2W3)\ + \ElWi]E{W2, W3)\ 

<4|p(Sfl,t)|+4|p(Sa,t)|. 

It then follows by induction that 

ka,t 







1 ka,, \ 


\ 


E 


exp 










[ i=l I 





i=l 



as — ► 00. 



□ 
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Proof of Lemmd^ A simple induction argument shows that if Zi, . . . z,„, Wi, . . . , Wm £ C are of 
modulus at most 1, then 

m 



n - n 



1=1 



Wi\. 



Applying the remark to | E [exp [i^Ua^)) \ and, for large enough a, 1 1 - f D | , we see 

ka.t ki 



Yl^[exp[inUa,i]] 

!=1 



2ka,t 



■I'D 



E 

ka,t 
ka.t 



E[exp [i^Ua.i]] - 



,,2 



■tT) 



2ka,t 



-tD 



l-^^E[t/y+0(E[<J]-[l- 

\kaMUl,l]-t'D\ + \ka.tO[^[Ull] 



2ka,t 



Now, it follows from the spectral representation of D that | fco^jE j j - fD| ^0. 

Expanding j results in b\ ^ terms of the form Za,iZa,iZa,k- Suppose i < j < k and 
let Di - j - i and D2- k - j. Let C{a) - Clog a, where C > 4/3. We separate the terms in the 
expansion of ^ into one group containing those terms whose indices satisfy Di < da) and 
D2 < C{a); and another group with those terms for which Di > C(a) or D2 > C{a). Suppose 
Za,iZa,jZa,k is in the first group. Then 

E[Za,iZa,jZa,k] < E^'^ [|Z„,,f | E^'^* \\Za.jf] E^'^ [\Za,kf] = E[|Z„,,f ] = 0(fl/logT'fl). 

Moreover, of the b^^ total terms, 0[ba,tC[a) 
contribution of these terms to | j O |e | [7^ j j 
negligible as a ^ 00. 

Suppose next that Za,iZajZa,k is in the second group and assume without loss of generality 
that Di > C{a). Observe that 

E[Za,iZa,jZa,k\ = corr(Za,,-, Z„jZ„,fc) E^'^ ,] (E[(Z^jZ^,fc)'] -E2 [Z„jZ„,fc])'" . 



of them fall into this first group. Thus the 
is at most of the order O (log^"'' a) which is 



Recall that E |z^ ^ J - O (log a). Moreover, by the Cauchy-Schwarz inequality, stationarity, and 
further direct computation of moments 

E\[Za,jZa,kf\ < e[z^,.] = 0[a^l\oi^a] andE^ [Z„,;Z„,fc] < E^ [z^ .] = Opog^a). 

Finally, using the mixing properties of the process con [Za,i,ZajZa,k) - P (C(fl)) - 0[a~'^) . 
Putting these three estimates together, we see that 

a-^E[Za,iZa.jZa,k] < O (fl-'^^^) iQgl/2-r 

Since, 0[b^ ^ - ba,tC{a)^] - 0(fo^ J - O^a^^log"^^ terms are in this second group, their 
contribution to |fca,fO|E |[/^ j] j | will be of order at most o\a'^^~'"\og~^'^~^^~^ a^, which is 



indeed negligible with our choice of C. The lemma now follows. 



□ 
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